A second-order space–time accurate scheme for nonlinear diffusion equation with general capacity term

A fully implicit finite difference (FIFD) scheme with second-order space–time accuracy is studied for a nonlinear diffusion equation with general capacity term. A new reasoning procedure is introduced to overcome difficulties caused by the nonlinearity of the capacity term and the diffusion operator in the theoretical analysis. The existence of the FIFD solution is investigated at first which plays an important role in the analysis. It is established by choosing a new test function to bound the solution and its temporal and spatial difference quotients in suitable norms in the fixed point arguments, which is different from the traditional way. Based on these bounds, other fundamental properties of the scheme are rigorously analyzed consequently. It shows that the scheme is uniquely solvable, unconditionally stable, and convergent with second-order space–time accuracy in L^∞(L²) and L^∞(H¹) norms. The theoretical analysis adapts to both one- and multidimensional problems, and can be extended to schemes with first-order time accuracy. Numerical tests are provided to verify the theoretical results and highlight the high accuracy of the second-order space–time accurate scheme. The reasoning techniques can be extended to a broad family of discrete schemes for nonlinear problems with capacity terms.     

Fast Evaluation of Linear Combinations of Caputo Fractional Derivatives and Its Applications to Multi-Term Time-Fractional Sub-Diffusion Equations

In the present work, linear combinations of Caputo fractional derivatives are fast evaluated based on the efficient sum-of-exponentials (SOE) approximation for kernels in Caputo fractional derivatives with an absolute error $\epsilon,$ which is a further work of the existing results in [13] (Commun. Comput. Phys., 21 (2017), pp. 650-678) and [16] (Commun. Comput. Phys., 22 (2017), pp. 1028-1048). Both the storage needs and computational amount are significantly reduced compared with the direct algorithm. Applications of the proposed fast algorithm are illustrated by solving a second-order multi-term time-fractional sub-diffusion problem. The unconditional stability and convergence of the fast difference scheme are proved. The CPU time is largely reduced while the accuracy is kept, especially for the cases of large temporal level, which is displayed by numerical experiments.     

On inexact alternating direction implicit iteration for continuous Sylvester equations

In this paper, we study the alternating direction implicit (ADI) iteration for solving the continuous Sylvester equation AX + XB = C , where the coefficient matrices A and B are assumed to be positive semi‐definite matrices (not necessarily Hermitian), and at least one of them to be positive definite. We first analyze the convergence of the ADI iteration for solving such a class of Sylvester equations, then derive an upper bound for the contraction factor of this ADI iteration. To reduce its computational complexity, we further propose an inexact variant of the ADI iteration, which employs some Krylov subspace methods as its inner iteration processes at each step of the outer ADI iteration. The convergence is also analyzed in detail. The numerical experiments are given to illustrate the effectiveness of both ADI and inexact ADI iterations.     

Accelerated projection-based forward-backward splitting algorithms for monotone inclusion problems

In this paper, based on inertial and Tseng''s ideas, we propose two projection-based algorithms to solve a monotone inclusion problem in infinite dimensional Hilbert spaces. Solution theorems of strong convergence are obtained under the certain conditions. Some numerical experiments are presented to illustrate that our algorithms are efficient than the existing results.     

A block-by-block method for the impulsive fractional ordinary differential equations

In this paper, a block-by-block numerical method is constructed for the impulsive fractional ordinary differential equations (IFODEs). Firstly, the stability and convergence analysis of the scheme are established. Secondly, the numerical solution which converges to the exact solution with order $3+\gamma$ for $0<\gamma<1$ is proved, where $\gamma$ is the order of the fractional derivative. Finally, a series of numerical examples are carried out to verify the correctness of the theoretical analysis.     

Two‐grid method for two‐dimensional nonlinear Schrödinger equation by finite element method

A conservative two‐grid finite element scheme is presented for the two‐dimensional nonlinear Schrödinger equation. One Newton iteration is applied on the fine grid to linearize the fully discrete problem using the coarse‐grid solution as the initial guess. Moreover, error estimates are conducted for the two‐grid method. It is shown that the coarse space can be extremely coarse, with no loss in the order of accuracy, and still achieve the asymptotically optimal approximation as long as the mesh sizes satisfy in the two‐grid method. The numerical results show that this method is very effective.     

A higher order difference method for singularly perturbed parabolic partial differential equations

This paper studies a higher order numerical method for the singularly perturbed parabolic convection-diffusion problems where the diffusion term is multiplied by a small perturbation parameter. In general, the solutions of these type of problems have a boundary layer. Here, we generate a spatial adaptive mesh based on the equidistribution of a positive monitor function. Implicit Euler method is used to discretize the time variable and an upwind scheme is considered in space direction. A higher order convergent solution with respect to space and time is obtained using the postprocessing based extrapolation approach. It is observed that the convergence is independent of perturbation parameter. This technique enhances the order of accuracy from first order uniform convergence to second order uniform convergence in space as well as in time. Comparative study with the existed meshes show the highly effective behavior of the present method.     

Forbidden sets and stability in some rational difference equations

In this paper, we solve open problem (5) submitted by Sedaghat in his paper, On third order rational difference equations with quadratic terms, J. Differ. Equ. Appl., 14(8) (2008), pp. 889–897. We also confirm conjecture (6) in the mentioned paper.     

On the optimal convergence factor of the accelerated parameterized inexact Uzawa method with three parameters for augmented systems

Under the assumption that all eigenvalues of the preconditioned Schur complement are real, we present an analytical proof for obtaining the optimal convergence factor of the real accelerated parameterized inexact Uzawa (APIU) method when P=A. It is proved that the optimal convergence factor is the same as that of the generalized successive overrelaxation method, which was published at the same time, and that it can be attained only at the unique optimum point of parameters, regardless of whether m>n or m=n. In addition, we generalize the APIU method and analyze the relationship between the APIU method and 10 additional Uzawa‐like methods.